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Abstract. The Poynting vector, energy density and energy velocity of light pulses 
propagating in anomalous dispersion medium (used in WKD-like experiments) are 
calculated. Results show that a negative energy density in the medium propagates 
along opposite of incident direction with such a velocity similar to the negative group 
velocity while the direction of the Poynting vector is positive. In other words, one 
might say that a positive energy density in the medium would propagate along the 
positive direction with a speed having approximately the absolute value of the group 
velocity. We further point out that neither energy velocity nor group velocity is a 
good concept to describe the propagation process of light pulse inside the medium in 
WKD experiment owing to the strong accumulation and dissipation effects. 



PACS numbers: 42.50.Gy, 42.25.Bs 



Wang, Kuzmich and Dogaiu claim that they have observed the superluminal light propagation by 
sending the probe laser pulses through the particularly prepared atomic caesium vapour and measuring 
the deviation of the refractive index of the medium from that of the vacuum and the arrival time of the 
pulses||l[ ^. Their main conclusion is that the light pulse propagates through the medium with group velocity 
about — c/310. Some comments on the WKD experiment have been made [3-6] in different ways. Recently, 
we analyse the negative group velocity and distortion of light pulse in the anomalous dispersion medium with 
two gain lines and show that although the shape of the light pulses seems to be almost preserved in a specific 
view, the distortion is still too large to use the concept of the group velocity to describe the propagation 
of the light in the mediumQ. In the present letter, we shall study the Poynting vector, energy density and 
energy velocity of light pulse in the anomalous dispersion medium. 

Without loss of generality, let the incident light pulse go to the medium at z = from vacuum in z < 0. 
The electric and magnetic fields E(r, and H(r, i) of the light pulse propagating in the medium in the 
z-direction can be expanded by complex Fourier series as 



E(z,t) = ( [ dtjE(z,tj)e-*('^+"<=)* +C.C.) 

2v27r \J ) 



and 



H(z,i) = ( / dcjH(z,w)e' 

2v27r \J 



(1) 
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where Wc is the angular frequency of the carrier wave, uj is the difference between the angular frequencies 
of the Fourier component and the carrier wave, E(z,a;) and H(z,a;) are the complex Fourier components 
of the electric and magnetic fields, respectively, with the angular frequency w + Wc, and c.c. stands for the 
complex conjugate. For monochromatic plane wave with angular frequency w + Wc, we have 



H(z, u) = VeRefc X E(z, w), (3) 
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where e{uj) is the dielectric permeabihty, and is the unit vector of wave. We have here chosen the magnetic 
conductivity ^ = 1 and then the refractive index is n{uj) = y/e{ij). Note that E, H and efe are perpendicular 
to each other and form a right-handed system. 

Substituting the above equations in the Poynting vector 

S = ^(ExH) (4) 

and taking the average over a period of carrier wave, we obtain 

S = ^e-^^MHo))- dwA{z,uj)e-'"' j (ko'n*{io')A*{z,u,')e''''' + c.c^ , (5) 

where fc(0) is the wave number of the carrier wave, and A{z,uj) = es ■ E(^;,a;)e~''^^°^^. Note that 

A{z, u) = A{0, u)) exp[i(fc(w) - k{0))z] (6) 
with k{uj) = {lu + LUc)n{ijj)/c and that 

k{ui) - k{0) = Kiw + ^K2u'^ + (7) 

where ki = (9'fc(u;)/9w')o with Z = 1,2, .... The subscript stands for taking value at zero-deviation from 
the carrier frequency. Thus, in order to obtain the Poynting vector, one has to calculate the integrals 

Ii = —z= / dujAiO,Lij) exp[—iLijt + iKiLoz + iK2(jj'^z/2 + ...] (8) 
v2n J-oo 

and 

1 f^°° 

I2 = —1= i dhjn{w)A{Q,w) ey.])[—iujt + iKiioz + in^ijJ^ z/2 + ...]. (9) 
v27r J-oo 

Now, consider the Guassian pulse, whose envelop at ^; = has the form 

A(0,t) = Aoexp[-a2(t-io)V2], (10) 

where a and to are two real constants, which determine the width and maximum position of the pulse 
respectively. Its Fourier transformation is 

A(0, t^) = ^ exp[-a;V(2a')]. (11) 



Keeping k,2 term and neglecting the higher order terms, one gets 

{t — tQ — KizY 



h — — = exp 



lK20fiz 



2 1 — iK20l Z 



(12) 



^2=5/1, (13) 

where 

g = n{uJs) + -- (14) 

^ 1 — lK2a z 

in which n"(0) = {d'^n(uj) / duj'^)^^^^ is the second derivative of refractive index with respect to w at zero- 
deviation from the carrier frequency, and 

.a^{t - to - mz) , , 

'^s = : 2 — -■ (15) 

1 — IK2OL z 
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Thus, the mean Poynting vector is 
5i cefe 



Stt |1 

The intensity of hght is given by 



° Re{g) exp 



-2Ini[fc(0)]z - a^Re 



1 — iK,2CX^Z 



(16) 



(17) 



where A is the area of the cross-section of the hght beam. It is evident that the intensity of hght / > as 
long as Re(5) > 0. When Re((7) ~ 1, the intensity of hght reduces to 



1 



|1 — iK2a z 



exp 



-a^Rc 



{t-tg- KizY 
1 — iK2Ct'^Z 



(18) 



which is the same as Eq.(31) in except the definition of ki here is equal to the one in [Q plus 1/c. 
In WKD experiment^ ||, cj^ = Ttt x 10*^MIIz, a 0.45MIIz and 



n{uj) 



1 + Anm 



1 



1 



LU — LUq + if] LU + LUo + if] 



1/2 



(19) 



where ujq = 2.77rMHz is a half of the angular frequency space of the two gain lines, (3 = 0.927rMHz is the 
angular frequency width of the gain lines, and m = O.SGttIIz is related to the gain coefficients. Neglecting 



m 



(20) 



and higher order terms, we obtain||l ki « -1.026 x 10"^ - 5.276 x lO^^''^^ s(cm)"\ K2 « 1.965 x lO"^^ _ 
2.415 x lO^^^i s^(cm)~^, and Re((7) = 1. Thus, the direction of the energy current in the medium is the 
same as the incident one of the pulse, and Eq.(p^ is a good approximate expression for the intensity of light. 

Now, let's turn to study the electromagnetic energy density and energy velocity in the medium. Using 
Maxwell's equations, the time average of divergence of Poynting vector over a period of the carrier wave can 
be written as 



V • S = — I E- 

47r 



'dt 



H 
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'dt 



dT> 1 9(H.H) 

47r i dt 2 dt 



Similar to E(z, t) and H(z, <), the induction D(z, t) can be expanded by complex Fourier series as 

1 



D(z,i) 



2V27r 



duje{Lo)'Ei{z, ijj)e 



+ c.c. 



Thus, 



'dt 



-zf{0)A{z,t) + f'iO) 



dA{z,t) 
dt 



2/ W 



+ ... ] e-*"'^*+*'=(")^ + C.C. 



(21) 



(22) 



(23) 



where /{lo) — (uj+ujc)e{uj), the prime stands for the derivative with respect to w, and A.{z, t) = £(2:, t) exp(iwci- 
ik{0)z) is the complex envelop of electric field strength E(z,t). The time average of the first term on the 
right-hand side of Eq. (Ell) in a period of carrier wave is 



9D 1 9[Re[/'(0)]A(z,t).A*(z,^)e-^i"'[''(")l-] , ^ ^ 
ar = 4 dt 



(24) 



where 



Q = ^ |2c^Jm(e(0))|A(z,t)|2 + 21m[/'(0)]lm (^A{z,t) 

d^A{z,ty 



Im 



A(z,i). /"(O) 



di' 



dA*{z,t) 
di 

I „-2Im[fc(0)]2 



(25) 
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Substituting Eq.(||) with ^) in (|2l|), we have 

dw 



dt 



V-S = -Q, 



with 



Stt 



iRe[/'(0)] A(z, t) . A* (2, i)e-2i-P(0)l- + H(z, i) • H(z, i)) 



(26) 



(27) 



It should be noted that the bars denoting the time average can be and hence have been moved inside the 
derivative operators on the left-hand side of Eq.(^). It should also be noted that in Eq.(^), A{z,t) is 
the complex envelop of the electric field, while H(z, t) is the real strength of magnetic field given by Eq.(||). 
Eq. ( p6|) is the differential conservation equation of the electromagnetic energy, w is the difference between the 
mean internal energy per unit volume with and without the electromagnetic field It may be interpreted 
as the electromagnetic energy density in the medium when the light pulse gets through. Q is the mean heat 
evolved per unit time and volume, representing the absorption (or dissipation) rate of energy per unit volume. 
The negative absorption (or dissipation) means that there exits gain (or accumulation) in the medium. One 
can define the velocity of transfer of energy density (i.e. energy velocity) as 



For the Guassian pulse Eq. ([lo|) , as an example 
_ Al Re[/'(0)] + 



167r \l-iK2a^z\ 



■ exp 



-2lm{k{Q))z - a^Re 



1 — iK2C? Z 



(28) 



(29) 



It can be checked that Eq.(26) holds for the Guassian pulse with Eqs.(|l^), ([l6| ) and ( ^9| ) if the terms in 
Eq.(^o|), the terms containing /'"(O) and higher order terms are neglected. In the approximation, 



/ d[(w + Wc)Re(n)/c 



Namely, the energy density transfers with the velocity approximately equal to group velocity. 
For WKD experiment. 



(30) 



Re[/'(0)] + « Re(n)Re 



duj 



-310 < 0, 



(31) 



so that w < 0. If one interpreted id as the electromagnetic energy density of the pulse in the medium, 
then the electromagnetic energy density in the medium would be negative! Besides, the direction of increase 
of energy density in the medium would be the opposite of that of |A(2;,t)|! The negative energy velocity 
Ve K Vg < would then imply that the negative energy density propagates in the opposite of the incident 
direction. 

To recover the positivity of energy density, one can replace tjj in Eq. 



I) by 



with 



Uq — const. > 0. 



(32) 



Uq is the zero-point energy density coming from the preparation of the medium so that the atomic population 
inverted, the part of electromagnetic energy of two beams for preparing the system being stored in caesium 
atoms. The probe pulse through the medium will cause that the number of inverted atomic population, so 
as the total energy density, decreases at some time. This is in accordance with the fact that the negative 
oscillator strength is assigned to obtain Eq.([l9|) ||^, |l^. Uq also sets the upper limit of the intensity of light. 
It reflects the possible saturation effect |^. The negative energy density propagating in the opposite of the 
incident direction is equivalent to the positive energy density propagating with |t;e| < c along the incident 
direction at any z in the medium. Mathematically, one may define the energy velocity in this case as 



S-efc 

Uo-U 



-Vg > 0. 



(33) 
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It implies that the propagation of positive energy density is subluminal at every point. 

It seems that there is a contradiction between the subluminal propagation of positive energy density 
and the peak of pulse arriving at the exit of the medium earlier than through the vacuum, even earlier than 
arriving at the entrance of the medium. However, the contradiction would exist only if Q could be negligible. 
In fact, Q in WKD experiment cannot be neglected. For the Gaussian pulse Eq.(|l0|), 

dw 

— = 21it,{los)w. (34) 
Im(a;s) depends on t and z. Its typical value for the pulse is about 10^s~^. While 

2c.Jm(.(0)) +Im[/"(0)]Re ( "\ ) + Im[/"(0)]Re 1 , (35) 

\l-iK.2a zj J Rc(/ (0)) + Igl 

where the second term in Eq.(|25|) has been neglected because it is much smaller than the other two terms. 
The first term in the braces in Eq.(|35|) gives a constant gain (or accumulation) rate, which is three orders 
of magnitude greater than the typical value of lOs- It magnifies the intensity of light by exp{— 2Im[fc(0)]2:} 
times. The second term, depending on z, and the third term, depending on both t and z, contribute the gain 
(or accumulation) and absorption (or dissipation) rates, respectively, and are about two orders of magnitude 
greater than the typical value of LOg- It is the latter two terms that are responsible for the advancement shift 
of the pulse. The second term shows that the gain rate increases as z increases. The third term shows that 
the absorption rate takes the minimum value at t = to + Kiz. Although the absorption rate is time-symmetric 
around t = to + kiz, the effect is very different because the intensity of light is very small when t < to + niz 
and becomes larger when t > to + kiz. In words, the gain (or accumulation) rate in WKD-like experiments 
are not constant as expected in 

Combining the study on the propagation of the wave front with a step- function in we can now give 
the following picture for the WKD experiment: The leading edge of the Guassian pulse entering the medium 
induces the gain (or accumulation) effect that makes the peak of the pulse to appear at the exit of the 
medium earlier than it through the vacuum. At later time when the peak of the pulse reaches the medium, 
the absorption (or dissipation) effect becomes more important, and so that a part of the pulse energy would 
be dissipated inside the medium. Thus the peak at the exit of the medium should not be the entrance 
peak. This picture is in agreement with Zhang's argument about gain effect pl , McDonald's energy loan 
argument 1^, as well as Sprangle-Pehano-Hafizi's differential gain argument ||^. In other words, the process 
of the pulse through the caesium vapour is not a simple propagation process of energy because there exists 
strong interaction between the light pulse and the medium. In this sense, the energy velocity as well as the 
group velocity are all not good concepts of describing the process of light pulse propagating through the 
medium in WKD experiment. 
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